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By means of neutron scattering measurements, we have investigated spin- wave excitation in 
a collinear four-sublattice (4SL) magnetic ground state of a triangular lattice antiferromagnet 
CuFeCb, which has been of recent interest as a strongly frustrated magnet, a spin-lattice coupled 
system and a multiferroic. To avoid mixing of spin- wave spectrum from magnetic domains 
having three different orientations reflecting trigonal symmetry of the crystal structure, we 
have applied uniaxial pressure on [110] direction of a single crystal CuFeCb- By elastic neutron 
scattering measurements, we have found that only 10 MPa of the uniaxial pressure results 
in almost 'single domain' state in the 4SL phase. We have thus performed inelastic neutron 
scattering measurements using the single domain sample, and have identified two distinct spin- 
wave branches. The dispersion relation of the upper spin- wave branch cannot be explained by 
the previous theoretical model [R. S. Fishman: J. Appl. Phys. 103 (2008) 07B109]. This implies 
the importance of the lattice degree of freedom in the spin- wave excitation in this system, 
because the previous calculation neglected the effect of the spin-driven lattice distortion in the 
4SL phase. We have also discussed relationship between the present results and the recently 
discovered " electromagnon" excitation. 
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1. Introduction 

From 1990s, a delafossite compound CuFe02 has been 
extensively studied as a frustrated spin system. 1 ' 2 ) Be- 
cause of triangular geometry of magnetic Fe 3+ ions and 
antiferromagnetic interactions between them, this sys- 
tem has vast degeneracy around the ground state, in 
which a collinear four-sublattice (4SL) antiferromagnetic 
order is realized. ^ Therefore, application of magnetic 
fields or substitution of a small amount of nonmagnetic 
ions for magnetic Fe 3+ ions induces a variety of mag- 
netic phases and unconventional magnetic phase tran- 
sitions. 3 8 ) Moreover, the first magnetic-field induced 
phase from the 4SL phase has been recently focused 
as a spin-driven ferroelectric phase, in which an elliptic 
screw-type magnetic ordering breaks inversion symme- 
try of the system and accounts for the charge polarized 
state. 9 n ) To understand these various magnetic phase 
transitions, it is essential to establish the spin Hamilto- 
nian of this system. One of the most effective way to 
determine the spin Hamiltonian is to measure spin- wave 
dispersion relation by inelastic neutron scattering mea- 
surements. Although several previous studies have re- 
ported inelastic neutron scattering measurements in the 
4SL phase, 12 14 ) there still remains some ambiguity in 
determination of the spin- wave dispersion relations from 
the measured spectrum. This is due to existence of three 
magnetic domains reflecting trigonal symmetry of the 
crystal structure. Figures 1(a) and 1(b) show the crys- 



tal structure of CuFeC^ and the 4SL magnetic structure, 
whose magnetic propagation wave vector is ( - , - , - ) , re- 
spectively. Owing to the threefold rotational symmetry 
about the c axis, there are three magnetic domains whose 

wave vectors of (|,-^,f) and (~h h f) are 

crystallographically equivalent to each other. Therefore, 
the magnetic excitation spectrum measured in the previ- 
ous studies are superpositions of the spin- wave spectrum 
from the three magnetic domains with different orienta- 
tions. On the other hand, recent synchrotron radiation 
x-ray diffraction studies have revealed that CuFe02 ex- 
hibits trigonal to monoclinic crystal structural transition 
associated with the paramagnetic to antiferromagnetic 
phase transition. 15 ' 16 ^ This structural transition also re- 
sults in three monoclinic structural domains, and there 
is a one-to-one correspondence between the monoclinic 
structural domains and the magnetic domains, as shown 
in Figs. l(c)-l(e). This implies the possibility that vol- 
ume fractions of the three magnetic domains can be con- 
trolled by application of pressure, which directly affects 
lattice degree of freedom. In the present study, we have 
demonstrated that application of uniaxial pressure on a 
single crystal CuFe02 results in almost 'single-domain' 
4SL phase. Using the single-domain sample, we have per- 
formed inelastic neutron scattering measurements, and 
have identified the spin- wave dispersion relations for the 
4SL phase. 
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Fig. 1. (Color online) (a) Crystal structure of CuFeC>2 with the 
hexagonal basis, (b) Magnetic structure of the 4SL phase with 
the monoclinic basis, [(c)- (e)] Relationships between the hexag- 
onal and the monoclinic bases in the (c) (110)-, (d) (120)- and 
(e) (210)-domains. 



2. Preliminary Details and Experiments 

CuFe02 exhibits two magnetically ordered phases in 
zero magnetic field. The high temperature phase is an 
incommensurate collinear magnetic phase referred to as 
partially disordered (PD) phase, 17 ) whose propagation 
wave vector is (g, g, §) with q = 0.202 ~ 0.210, and the 
low temperature phase is the 4SL phase. The transition 
temperatures for the two magnetic phases are Tni = 14 
K and Tn2 ~ H K, respectively. Below Tni, the crys- 
tal structural transition occurs and the three magnetic 
domains are formed. In this paper, we refer to the three 



domains with the wave vectors of (|, |, §), (\->— ^? f ) 
and (-|,|,f) as ( 110 K (120)- and (210)-domains, re- 
spectively. Relationships between the monoclinic and the 
hexagonal bases in each domain are shown in Figs. 1(c)- 
1(e). To distinguish between the two bases, the subscript 
"m" has been added to the monoclinic notation when 
referring to the basis, wave vectors and reciprocal lattice 
indices. The previous x-ray diffraction studies have re- 
vealed that with decreasing temperature, a m contracts, 
and on the contrary, b m elongates in each domain. 15,16 ) 
We thus anticipated that uniaxial pressure applied per- 
pendicular to the c axis favors the magnetic domains 
whose a m axis lies along the pressure, and, on the con- 
trary, suppresses those whose b m axis lies along the pres- 
sure. Therefore, we cut a single crystal of CuFe02 grown 
by the floating zone technique 18 ) into thin plate shape 
(~ 6x10x2 mm 3 ) with the widest surface normal to the 
[110] direction, on which we applied uniaxial pressure so 
as to maximize the volume fraction of the (llO)-domain. 
The sample was set in a uniaxial pressure cell devel- 
oped by Aso et al. 19 ^ We applied, at room temperature, 
10 MPa of uniaxial pressure on the sample. CuBe disk 
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Fig. 2. (Color Online) (a) Comparison between the observed and 
the calculated magnetic structure factors. Open circles, triangles 
and squares denote data for the (110)-, (120)- and (210)-domains, 
respectively. Inset shows schematic drawing of direction of the 
uniaxial pressure (open blue arrows) and orientations of the three 
magnetic domains. Filled red, green and purple arrows denote di- 
rections of c-plane-projections of the magnetic modulation wave 
vectors; the sizes of the arrows qualitatively show the volume 
fractions of the three domains, (b) Temperature variation of in- 
tegrated intensities of magnetic reflections corresponding to the 
PD and the 4SL magnetic orderings in the (llO)-domain. 



springs were set in the pressure cell to keep the applied 
pressure. In order to evaluate the volume fractions of the 
three magnetic domains, we performed four-circle neu- 
tron diffractions measurement using FONDER installed 
at JRR-3 in Japan Atomic Energy Agency (JAEA). The 
incident neutron beam with wavelength 1.240 A was ob- 
tained by a Ge(311) monochromator. The sample with 
the pressure cell was mounted on a closed-cycle He-gas 
refrigerator, and was cooled down to 2.5 K. 

After the four-circle neutron diffraction measurement, 
we have performed inelastic neutron scattering measure- 
ment using the identical sample in the pressure cell, 
which has kept the 10 MPa of uniaxial pressure through- 
out the two neutron scattering experiments. For this 
measurement, we used a cold neutron triple axis spec- 
trometer HER(Cl-l) installed at JRR-3 in JAEA. The 
wave number of the scattered neutrons was fixed to 
kf — 1.3246 A -1 . A horizontal focusing analyzer was em- 
ployed. Energy resolution at elastic position is 0.13 meV 
(full width at half maximum). The higher-order contam- 
inations were removed by a cooled Be-filter and a PG fil- 
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Fig. 3. (Color Online) (a) Typical profiles of constant-Q scans in the present inelastic neutron scattering measurement. Open and filled 
allows denote peak positions corresponding to the HE and the LE branches, (b) Reciprocal lattice map of the (H, H, L) scattering plane. 
Open squares and a circle denote positions of nuclear and magnetic Bragg reflections, respectively. Filled squares denote positions of 
magnetic zone center, [(c)-(g)] Spin-wave spectrum along (c) (H , H , §), (d) (0,0, L), (e) (±, |,L), (f) (±, \,L) and (g) (0.29, 0.29, L) 
lines. Circles and squares represent peak positions in constant-Q scans, and vertical bars represent full width at half maximum of the 
peaks. Dashed line in (c) shows the asymmetric branch observed in the previous studies 13 ' 14 ) (taken from Ref. 14). 



ter placed in front of the sample and a room-temperature 
Be-filter placed in front of the analyzer. The sample with 
the pressure cell was mounted in a pumped 4 He cryostat 
with (i7, H, L) scattering plane, on which the magnetic 
modulation wave vector of the (llO)-domain lies, and was 
cooled down to 1.4 K. 

3. Results and Discussions 

3.1 Four- circle neutron diffraction measurements 

In the present four-circle neutron diffraction measure- 
ment, we measured intensities of 32 magnetic Bragg re- 
flections, of which 15, 7 and 10 magnetic reflections be- 
long to the (110)-, (120)- and (210)-domains, respec- 
tively. The effect of neutron absorption was corrected by 
the DABEX software. We have performed least-square 
fitting analysis, in which we used the established mag- 
netic structural parameters of the 4SL phase presented 
in Ref. 1 and refined the volume fractions of the (110)-, 
(120)- and (210)-domains, V(no), V(i20) an d V^io)- The 
comparison between the observed magnetic structure fac- 
tor | F b s | and the calculated values |-F C ai| is shown in 
Fig. 2(a). The reliability factor was obtained to be 7.2 
%, indicating that all of the observed data is well ex- 
plained by the 4SL magnetic structure. The volume frac- 
tions of the three magnetic domains are determined to be 
V(iio) • V(i2o) • V(2io) = 0.973 : 0.012 : 0.015. This shows 
that the uniaxial pressure applied onto the [110] surface 
successfully produces almost 'single-domain' 4SL state in 
which the (llO)-domain dominates over the (120)- and 
the (210)-domains. 

Figure 2(b) shows temperature variation of integrated 
intensities of magnetic reflections corresponding to the 



PD and the 4SL magnetic orderings in the (llO)-domain, 
which were measured on cooling. We found that the tran- 
sition temperatures and temperature variation of the or- 
der parameters are almost the same as those in zero pres- 
sure. 16 ) This indicates that the application of 10 MPa of 
the uniaxial pressure does not largely affect the magnetic 
interactions in this system, while it does the volume frac- 
tions of the three domains. 

3.2 Inelastic neutron scattering measurements 

Figures 3(a) and 3(c) show typical profiles of constant 
wave vector (Q) scans and magnetic excitation spectrum 
along (#, i^, |) line, respectively. These data have re- 
vealed that in the 4SL phase, there are two different 
spin- wave branches. In this letter, we refer to the upper 
and the lower branches as high-energy (HE) and low- 
energy (LE) branches, respectively. We found that both 
of the branches are symmetric with respect to H = ~. 
On the other hand, in the previous studies, 13, 14 ) another 
asymmetric branch has been observed in the (H,H, |) 
line; the dispersion relation of the asymmetric branch is 
drawn by a dashed line in Fig. 3(c). Ye et al. have sug- 
gested that the asymmetric branch belongs to the mag- 
netic domains whose propagation wave vectors point out 
of the (H, H, L) scattering plane. 14 ) Actually, the present 
measurements have demonstrated that the asymmetric 
branch was not observed in the 'single-domain' 4SL state, 
confirming the previous suggestion. 

It should be mentioned that recent theoretical work by 
Fishman has pointed out that the spin-wave dispersion 
along the (H, H, |) line consists of two branches. 20 ) How- 
ever, in the previous inelastic neutron scattering studies 
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on this system, only the LE branch has been identified 
to belong to the (llO)-domain. 14 ' 21 ) This is because the 
energies of the HE branch are close to those of the asym- 
metric branch, and therefore, it is quite difficult to distin- 
guish between them. Removing the asymmetric branch 
by the applied uniaxial pressure, we have, for the first 
time, identified the dispersion relation of the HE branch. 
This suggests that the magnetic domain control by ap- 
plied uniaxial pressure is useful to determine spin-wave 
dispersion relation in spin-lattice coupled complex mag- 
nets with multi-domain structures. 

The measured dispersion relation of the HE branch has 
two energy dips at H ~ 0.33 and ~ 0.17. This does not 
agree with the previous theoretical result by Fishman, 20 ) 
indicating that the magnetic interaction parameters pre- 
sented in Refs. 14 and 20 are not sufficient to explain the 
spin- wave dispersion relation in this system. We have also 
measured the spin-wave dispersion along the c* axis at 
H = K = 0.25 and 0.29, as shown in Figs. 3(f) and 3(g). 
These results have revealed that the HE branch is less 
dispersive as compared to the LE branch. 

We now discuss the energies of the spin waves at 
the magnetic zone center, specifically, (0,0,0), (0,0,3), 
(i> i> §)> an d so on 7 because these energies are 

relevant to optical responses including ESR modes. As 
for the LE branch, the previous work by Ye et al. 14 ^ have 
clearly showed that there are two peaks at 1.1 and 1.4 
meV in a constant-Q scan at (0, 0, 3) [see Fig. 2(b) in Ref. 
14]. Although in the present measurements, we could not 
identify the two-peak structure in the constant-Q scan at 
(0,0,3) due to lack of intensity, the asymmetric peak pro- 
file in the constant-Q scan at (|, |, |) suggests that the 
LE branch is split into two branches at the magnetic 
zone center [see Fig. 3(a)]. These energies show good 
agreements with the antiferromagnetic resonance modes 
around 270 and 340 GHz observed in the previous ESR 
measurements on CuFe02- 22 ' 23 ^ 

As for the HE branch, we have measured the spin- 
wave dispersion along the (|, ^, L) line as shown in Fig. 
3(e), and have found that the HE branch has an energy 
of 2.47 meV at (^,^,0). This energy agrees with the 
ESR signal around 600 GHz observed in the previous 
measurement. 23 ) However, in Ref. 23, Kimura et al. have 
argued that their conventional spin- wave calculation can- 
not reproduce the ESR signal around 600 GHz. On the 
other hand, quite recently, Seki et al. have performed 
terahertz time-domain spectroscopy in the 4SL phase. 
They have observed an 'electromagnon' (electric-field- 
active magnon) excitation around 2.3 meV, 24 ) which 
agrees with the energy of the HE branch at (^, |, | ). 
These results suggest that the HE branch has novel spin 
dynamics coupled with lattice or charge degree of free- 
dom. It should be also mentioned that the reason for the 
finite difference between the energy of the HE branch at 
(|, |,0) (2.47 meV) and that at (±, \, f ) (2.30 meV) is 
not clear at this moment. 

4. Conclusion 

In conclusion, we have investigated spin- wave excita- 
tion in the 4SL magnetic ground state of CuFe02- To 



avoid mixing of the spin-wave spectrum from the three 
magnetic domains, we applied 10 MPa of uniaxial pres- 
sure onto the [110] surface of the single crystal CuFe02- 
The present four-circle neutron diffraction measurements 
have demonstrated that the application of the pressure 
successfully produces 'single-domain' 4SL phase. Using 
the identical sample, we have performed the inelastic 
neutron scattering measurements. As a result, we have 
successfully identified the two (LE and HE) spin-wave 
branches, both of which are symmetric with respect to 
(l> \ -> §)• It should be emphasized that the dispersion 
relation of the HE branch does not agree with the pre- 
vious theoretical work, 20 ) which is based on the assump- 
tion that the spin-driven crystal lattice distortion 15 ' 16 ) 
hardly affects the spin-wave dispersion relations in the 
4SL phase. We have thus suggested that the spin-lattice 
coupling is essential to understand the spin-wave excita- 
tion in this system. The asymmetric branch observed in 
the previous works 12 14 ) was not observed in the single 
domain 4SL state. This confirms the previous suggestion 
that the asymmetric branch arises from the magnetic do- 
mains with the magnetic modulation wave vectors point- 
ing out of scattering plane. 14 ) We have also discussed the 
relationship between the energies of the spin waves at the 
magnetic zone center and the results of the optical mea- 
surements on this system, 22 24 ) suggesting that the HE 
branch corresponds to the "electromagnon" excitation, 
in which dynamical coupling between spin and dielectric 
polarization is realized. 

Acknowledgment 

We are grateful to Shojiro Kimura for fruitful discus- 
sions. The neutron scattering measurements at JRR-3 
were carried out along the proposals No. 10783B and 
No. 10586, and were partly supported by ISSP of the 
University of Tokyo. This work was partly supported by 
Grants-in-Aid for Scientific Research gYoung Scientists 
(B), Grants No. 20740209h. We thank M. Yokoyama for 
technical support in the neutron scattering experiment 
with HER(Cl-l). The images of the crystal and mag- 
netic structures in this letter were depicted using soft- 
ware VESTA. 25 ) 



1) S. Mitsuda, H. Yoshizawa, N. Yaguchi and M. Mekata: J. Phys. 
Soc. Jpn: 60 (1991) 1885 

2) M. Mekata, N. Yaguchi, T. Takagi, T. Sugino, S. Mitsuda, H. 
Yoshizawa, N. Hosoito and T. Shinjo: J. Phys. Soc. Jpn: 62 
(1993) 4474 

3) S. Mitsuda, M. Mase, K. Prokes, H. Kitazawa, and H. A. Katori: 
J. Phys. Soc. Jpn: 69 (2000) 3513 

4) N. Terada, S. Mitsuda, T. Fujii, K. Soejima, I. Doi, H. A. Katori 
and Y. Noda: J. Phys. Soc. Jpn. 74 (2005) 2604. 

5) N. Terada, Y. Narumi, K. Katsumata, T. Yamamoto, U. Staub, 
K. Kindo, M. Hagiwara, Y. Tanaka, A. Kikkawa, H. Toyokawa, 
T. Fukui, R. Kanmuri,T. Ishikawa and H. Kitamura: Phs. Rev. 
B 75 (2007) 224411. 

6) S. Seki, Y. Yamasaki, Y. Shiomi, S. Iguchi, Y. Onose, and Y. 
Tokura: Phys. Rev. B 75 (2007) 100403(R). 

7) N. Terada, T. Nakajima, S. Mitsuda, H. Kitazawa, K. Kaneko, 
and N. Metoki: Phys. Rev. B 78 (2008) 014101. 

8) E. Pachoud, C. Martina, B. Kundys, C. Simona, and A. Maig- 
nana: J. Solid State Chem. 183 (2010) 344. 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 5 



9) T. Kimura, J. C. Lashley and A. P. Ramirez: Phs. Rev. B 73, 

(2006) 220401(R). 

10) T. Nakajima, S. Mitsuda, S. Kanetsuki, K. Prokes, A. 
Podlesnyak, H. Kimura, and Y. Noda: J. Phys. Soc. Jpn. 76 

(2007) 043709. 

11) T. Nakajima, S. Mitsuda, K. Takahashi, M. Yamano, K. Ma- 
suda, H. Yamazaki, K. Prokes, K. Kiefer, S. Gerischer, N. Ter- 
ada, H. Kitazawa, M. Matsuda, K. Kakurai, H. Kimura, Y. 
Noda, M. Soda, M. Matsuura, and K. Hirota: Phys. Rev. B 
79 (2009) 214423. 

12) N. Terada, S. Mitsuda, Y. Oohara, H. Yoshizawa, and H. Takei: 
J. Magn. Magn. Mater. 272-276S (2004) E997. 

13) N. Terada, S. Mitsuda, T. Fujii, and D. Petitgrand: J. Phys.: 
Condens. Matter 19 (2007) 145241. 

14) F. Ye, J. A. Fernandez-Baca, R. S. Fishman, Y. Ren, H. J. Kang, 
Y. Qiu, and T. Kimura: Phys. Rev. Lett. 99 (2007) 157201. 

15) N. Terada, S. Mitsuda, H. Ohsumi, and K. Tajima: J. Phys. 
Soc. Jpn. 75 (2006) 023602. 

16) F. Ye, Y. Ren, Q. Huang, J. A. Fernandez-Baca, P. Dai, J. W. 
Lynn, and T. Kimura: Phs. Rev. B 73 (2006) 220404(R). 

17) S. Mitsuda, N. Kasahara, T. Uno, and M. Mase: J. Phys. Soc. 



Jpn. 67 (1998) 4026. 

18) T. R. Zhao, M. Hasegawa and H. Takei, J. Cryst. Growth 166 
(1996) 408. 

19) N. Aso, Y. Uwatoko, H. Kimura, Y. Noda, Y. Yoshida, S. Ikeda, 
and S. Katano: J. Phys.: Condens. Matter 17 (2005) S3025. 

20) R. S. Fishman: J. Appl. Phys. 103 (2008) 07B109. 

21) In Ref. 13, Terada et al. have presented the spin- wave branch 
which is also referred to as "HE-branch". However, the "HE- 
branch" in Ref. 13 is not symmetric with respect to the mag- 
netic zone center, and is entirely different from the 'true' HE- 
branch presented in this paper. This misidentification of the 
spin-wave dispersion relation is due to the contamination of 
the scattering from the other domains. 

22) T. Fukuda, H. Nojiri, M. Motokawa, T. Asano, M. Mekata, 
and Y. Ajiro: Physica B 246-247 (1998) 569. 

23) S. Kimura, N. Nishihagi, H. Yamaguchi, T. Kashiwagi, M. Hagi- 
wara, N. Terada, and K. Kindo: J. Low Temp. Phys. 159 (2010) 
80. 

24) S. Seki, N. Kida, S. Kumakura, R. Shimano, and Y. Tokura: 
Phys. Rev. Lett. 105 (2010) 097207. 

25) K. Momma and F. Izumi: J. Appl. Crystallogr. 41 (2008) 653. 



